Effect of self-affine fractal characteristics of surfaces on wetting 



S. Sarkar a , S. Patra a *, N. Gayathri 5 t, g. Banerjee c * 
a Department of Electronics and Communication Engineering, 
Haldia Insitute of Technology - Haldia West Bengal 721 657 India 
b Material Science Division, 
Indira Gandhi Center for Atomic Research Kalpakkam 60S 102 India 
c Surface Physics Division, 
Saha Institute of Nuclear Physics 1/AF Bidhannagar Kolkata 700 064 India 

The relation between the contact angle of a liquid drop and the morphological parameters of 
self-affine solid surfaces have been investigated. We show experimentally that the wetting property 
of a solid surface crucially depends on the surface morphological parameters such as: (1) root mean 
square (rms) roughness a, (2) in-plane roughness correlation length £ and (3) roughness exponent 
a of the self-affine surface. We have shown that the contact angle monotonically decreases with the 
increase in the rms local surface slope p (oc cr/^ a ) for the cases where the liquid wets the crevices of 
the surface upon contact. We have shown that the same solid surface can be made hydrophobic or 
hydrophilic by merely tuning these self-affine surface morphological parameters. 
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Surface morphology plays a very crucial role in deter- 
mining the wettability of a surface^. Knowledge about 
wetting in the sub-microscopic level will have immense 
impact in narrow channels fluid-flow technology^ and 
in biological applications^. Wetting on a smooth and flat 
solid substrate could be well described by Young's^ equa- 
tion and on rough surfaces of solid could be well described 
by Wenzel'*^ or Cassie-Baxter's 8 formalism. Wettabil- 
ity of a surface is mainly determined by three energies: 
the solid/vapor (SV) interfacial energy ^sv, solid/liquid 
(SL) interfacial energy jsl and liquid/vapor (LV) inter- 
facial energy jlv- Young pointed out the relationship 
between the contact angle and the above parameters. 
The Young's equation is expressed as: cosOy = (jsv — 
1sl)/ilv or Isl = Isv — JlvcosOy where Y is the con- 
tact angle the liquid droplet makes with the smooth and 
flat solid surface. When j S v > Isl then 0° < Y < 90° 
and for jsl > 7sv then 90° < Y < 180°. Introduction 
of roughness on solid surfaces, increases its interfacial 
energy with liquid and vapor. Wenzel 7 proposed that 
the modified contact angle due to roughness can be ex- 
pressed as cosOw = r(jsv — 1sl)/jlv = rcosOy where 
"r" is the ratio between the true area A and the pro- 
jected area Aq (r = A/Aq) and is always greater than 1. 
Cassie and Baxter^ further extended this by considering 
the solid/liquid (SL) interface to have air pockets in the 
crevices below the liquid drop due to high solid asperi- 
ties (features like -"Lotus effect"), where some fraction of 
the solid f$ is in contact with the liquid drop and some 
fraction fy under the drop has air or vapor pockets. By 
including the area ratio factor "r", Cassie and Baxter 
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showed the dependence of the modified contact angle as 
a function of f$ i.e., cosOcb = rfscosOy + /s — 1. Here, 
as fs tends to zero, the contact angle approaches 180° 
and as fs tends to 1, the expression tends to the Wenzel's 
equation. Recently we have shown a simple derivation of 
Young, Wenzel and Cassie-Baxter equation starting from 
the work of adhesion 9 . 

The area ratio factor " r" depends mainly on three mor- 
phological parameters (1) root mean square (rms) rough- 
ness a, (2) in-plane roughness correlation length £ and (3) 
roughness exponent a. If these parameters are isotropic, 
then the morphology is termed as self-affine, has a fractal 
dimension and is generally found to be valid for various 
growth process such as vapour deposition of metal films, 
fractured surfaces and many other natural surfaces. The 
above mentioned formalisms that explain wetting does 
not consider the surface morphological parameters such 
as growth roughness exponent a and lateral correlation 
length £. The above formalisms also do not tell us how 
to modify a surface from hydrophobic to hydrophilic or 
vice versa by simply manipulating parameters a and £. 
To our knowledge no experimental work has been car- 
ried out to investigate the relationships between wetting 
and surface morphological parameters a and £ but some 
theoretical attempt has been made 10 . In this paper we 
will present some experimental results showing very in- 
teresting behavior of wetting depending strongly on the 
value of these three morphological parameters and per se 
we can tune the wettability of any solid surface without 
changing the chemical composition. 

We have used the following definitions for extract- 
ing these parameters from the data: For a given sur- 
face, if h(x,y) is the height at position (x,y), then 
one can obtain the root mean square deviation of the 
height which is commonly termed as rms roughness 
a = (< [h(x,y)— < h(x,y) >] 2 >) 1//2 where the notation 
< ... > means the spatial statistical average. But, more 



quantitative information of the surface morphology can 
be extracted from the height-height correlation function 
H(R) defined as H(R) =< [h(x, y) - h(x' y')} 2 >= 2a 2 - 
2 < h(x, y)h(x', y') > where R = y/[(x - x') 2 + (y - y') 2 } 
and < h(x,y)h(x' ,y') >= C(R) is the autocorrelation 
function (ACF). For self-affine fractal random surface 
model 11 the ACF C(R) = a 2 exp[- (R / £) 2a ] where £ and 
a are the parameters as defined earlier. In this formal- 
ism a describes the fractal characteristcs of the random 
surface and is related to Hausdorff dimension Df by 
Df=d — a,d being the embedded dimension. When 
a = 1, C(R) is the Gaussian function, representing ran- 
dom surfaces of Gaussian correlation. 

Now let us introduce a ratio of rms roughness ampli- 
tude cr (out-of-plane roughness) and in-plane roughness 
correlation length £ i.e., cr/£. This ratio is the measure 
of slope of the long wavelength undulations on the sur- 
face and is also called the rms long wavelength roughness 
raticP^. Since we know that a is the measure of the de- 
gree of sharp local surface irregularity then one can also 
define a rms local surface slope p oc cr/£ a which repre- 
sents the measure of the slope of the short wavelength 
undulationf^ (i.e., the slope of the sharp local irregu- 
larity). Here we shall assume that the liquid wets the 
crevices on the surface upon contact and this is possible 
only for the weak roughness surfaces^ i.e., for a > 0.5 
and cr/£ <C 1. We would like to see in this experimental 
study whether the contact angle has any dependence on 
this parameter p which is a function of cr, £ and a. 

To obtain surfaces of different morphology, we de- 
posited ZnO thin films layer- by- layer on glass substrates 
by combining a simple sol-gel and spin-coating tech- 
nique. For this investigation we prepared six samples (a) 
2-layered, (b) 4-layered, (c) 8-layered sample annealed 
at 450° C and (d to f) 8-layered samples annealed at 
500°C, 550°C and 600°C respectively. The details of 
the sample preparation and characterisation are reported 
elsewhere^. We have used Atomic Force Microscope 
(AFM) in the contact mode to map the morphology of 
the sample surface. 

The wetting properties of these surfaces were inves- 
tigated using a home made sessile drop microscope. A 
drop of double distilled deionized water was placed on 
the sample surface using a micro syringe. The image of 
the drop formed was analysed using ImageJ software^. 
The experiment was repeated at different regions over 
the sample surface. In fig. 1 we show a typical captured 
image of the 5 pi water droplet along with the contact 
angle measurement for all the six samples. Two impor- 
tant features can be immediately observed. (1) With the 
increase in the number of layers for the 450° C samples 
the contact angle progressively increases and the samples 
become more hydrophobic and (2) with the increase in 
annealing temperature the contact angle decreases and 
the samples become hydrophilic. To understand this be- 
haviour we imaged the surface morphology of our samples 
using an AFM. 

In fig. 2 we show the AFM images of all the six sam- 
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FIG. 1: (Colour online) CCD captured images of the 5 pi wa- 
ter droplets along with the contact angle measurement indi- 
cated by green lines for all the six samples. The double-sided 
arrow indicates the scale. 




FIG. 2: (color online) Atomic Force Microscope (AFM) im- 
ages of all the six samples, a, a and £ are also indicated. 
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FIG. 3: (a,b) Height-height correlation function H(R) versus 
R, the value of a can be obtained from the fit. (c,d) we show 
the ACF C(R) as a function of R and the solid line shows the 
fit to obtain the parameter £. The arrows indicate the cross 
over of the roughness exponent around R ~ 0.02 /im 



pies. On direct observation of the AFM images of the 
450°C annealed samples (samples (a) to (c), 2-layer, 4- 
layer and 8-layer respectively), we find that the surface 
morphology/asperities become more spiky with increas- 
ing number of layers, but the surface remains more or less 
flat. Whereas, on increasing the annealing temperature 
the surface does not remain flat and large surface undu- 
lations are formed (see the z-axis scale in fig. 2). The 
top surface morphology depends on the growth condi- 
tion, and playing with these growth conditions we could 
obtain a set of six different surfaces having different mor- 
phological parameters. 

To obtain the quantitative information of the surface 
morphology we have extracted all the three parameters 
i.e., (1) <j, (2) £ and (3) a using the formalism de- 
scribed earlier. In figs. 3(a,b) we have plotted the 
height-height correlation function H(R) versus R and 
since H(R) oc R 2a when R « £ the value of a can 
be obtained from the fit. We also observe that H(R) 
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FIG. 4: (Color online) (a) 3D plot of contact angle versus in- 
plane correlation length £ and the roughness exponent a. (b) 
contact angle versus rms surface roughness The error in the 
contact angle is represented by the size of the data symbols 
(dots). 



tends to 2<j 2 when R » £ as expected for the self-affine 
surface morphology. It is interesting to note that for 
sample (c) [450° C, 8-layers] we observe a crossover of the 
roughness exponent around R ~ 0.02 \im indicating that 
the surface morphology has an hierarchical structure. It 
has been predicted theoretically that this sort of heirar- 
chial structure will exhibit dewetting property of the solid 
surface^ 4 . We could also obtain the value of the lateral 
correlation length £ using the above mentioned formal- 
ism and by using the value of a obtained from H(R). In 
figs. 3(c,d) we show the ACF C(R) as a function of R 
and the solid line shows the fit from which the parameter 
£ has been obtained. We observe that for sample (c) the 
surface appears spiky (i.e., sharp local irregularities) and 
flat, leading to low a and £ similar to the morphology 
proposed by Cassie-Baxter. This leads to a high con- 
tact angle of the liquid drop. We see that the sample 
(c) does not obey the weak surface roughness criteria as 
discussed above and hence this sample (c) will not scale 
with the rms local surface slope p. For sample (f) even 
though the roughness exponent a is low indicating large 
local slope of the sharp irregularities, the rms roughness 
amplitude a and the in-plane lateral correlation length £ 
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is larger indicating large undulations on the surface and 
the contact angle of the water drop is observed to be the 
lowest. This aspect will be addressed further below. In 
fig. 4(a) we show a 3D plot of the contact angle versus 
in-plane correlation length £ and the roughness exponent 
a. We clearly observe that the contact angle is a non- 
monotonic function of a and £. In fig. 4(b) we observe 
another interesting feature i.e., we observe that as the 
rms surface roughness a increases the contact angle de- 
creases. Now we would like to investigate whether there 
exists any closed functional form exits between the con- 
tact angle and all these surface morphological parameters 
a, a and £. 

From the above experimental observation we can infer 
that one can tune the contact angle by controlling or tun- 
ing the surface morphological parameters. Large contact 
angle (hydrophobic) for sample (c) could be explained 
by invoking Cassie-Baxter formalism, but how do we ex- 
plain the hydrophilic nature obtained for sample (f) hav- 
ing moderately low roughness exponent a but relatively 
large lateral correlation length £ because of prominent 
undulations (mounds/peaks and troughs/valleys) on the 
sample surface (see fig. 2(f)) leading to topography driven 
wetting? Seperation between the peaks and valleys (the 
undulations) is the lateral correlation length £. For sam- 
ple (f) the asperities are riding over the undulations. This 
can be observed directly from the AFM image and was 
also observed from the fast Fourier transform of the im- 
age. If undulations are introduced on the surface then the 
liquid will follow the undulations and will flow/spill-over 
into the trough/ valleys (crevices) unlike in the Cassie- 
Baxter case where the liquid cannot break the liquid film 
and fill the air-pockets between the dense asperities on 
a flat substrate. The spill-over of liquid into the adja- 
cent trough/valleys can happen because large fraction 
of these valleys will be connected because of large sur- 
face rms roughness (as can be seen from the z-scale of 
fig. 2(f)). The penetration of the liquid to this adjacent 
troughs/ valleys is very similar to the model (penetration 
model) proposed by Ishino et al 15 . Large value of a can 
only give rise to the formation of large troughs/ valleys 
such that the liquid can spill-over into these. For the 
case of sample (f) event hough the value of a is low, we 
obtain from our analysis a large value of the slope of the 
long wavelength undulation on the surface (cr/£), thus 
preventing the surface to go to the Cassie-Baxter limit. 
These large undulations help in transforming the hydrob- 
hobic surface to a hydrophilic surface. Thus, our experi- 
mental results indicate that all the three surface morpho- 
logical parameters cr, £ and a are important to determine 



the wetting property of the solid surface. 

Next, we tried to see if the contact angle has any de- 
pendence on these three surface morphological parame- 
ters a, £ and a. In fig. 5 we have plotted the contact 
angle versus rms local surface slope p. We observe that 
the contact angle has a monotonic dependence on p ie., 
the contact angle decreases with increasing value of p. 
The sample (c) is shown as an isolated point since it is 
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FIG. 5: Contact angle verses p (oc cr/^ a ) the rms local surface 
slope. The surface of Sample C is in the Cassie-Baxter limit. 



a true Cassie-Baxter surface. From fig. 5 it appears that 
the dependence of contact angle on p is valid only for the 
cases where the liquid wets the crevices of the surface 
upon contact. 

To summarize, in this report we have characterized the 
surface morphology of sol-gel grown ZnO multilayer films 
using atomic force microscopy and obtained the surface 
morphological parameters. These parameters were corre- 
lated to the wetting property (contact angle) of the liquid 
over the sample surface. We found that the contact angle 
monotonically depends on the rms local surface slope p 
for a self-affine fractal surface when the liquid wets the 
crevices of the sample surface, p is a function of all the 
three surface morphological parameters p(<r, £,a). We 
could infer from our experimental results that we can 
tune the wetting property (of self-affine fractal surfaces) 
by merely changing the morphological parameters of the 
surface. 
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